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[0001] The benefit of the Filing date of provisional U.S. application Serial Number 

60/554,456, filed 19 March 2004, is claimed under 35 U.S.C. § 119(e) in the United States, 
and is claimed under applicable treaties and conventions in all countries. 

[0002] The development of this invention was partially funded by grants 

R03EY0 12887 and P30EY02377 from the National Eye Institute, National Institutes of 
Health, Bethesda, Maryland; and from a Space Product Development grant from the National 
Aeronautics and Space Administration. The Government has certain rights in this invention. 

TECHNICAL FIELD 

[0003] This invention pertains to a method to measure relative changes in blood 

oxygen saturation using hemoglobin spectral curves generated from reflected light from in 
vivo blood vessels, e.g., retinal macro- and micro-circulation. 

BACKGROUND ART 

[0004] The oxygen supply of the retina is provided by both the choroidal and retinal 

circulation. Because of the high oxygen needs of the retina, any alteratipn in circulation such 
as seen in diabetic retinopathy, hypertension, sickle cell anemia, and vascular diseases can 
result in impairment. Pathological conditions in the retina and optic nerve head (ONH) can 
cause vision loss and blindness. Both structures have a high demand for oxygen, and loss of 
the normal oxygen supply through vascular insufficiency is believed to play an important role 
in retinal and ONH pathology. See, G. A. Cioffl et al, "Optic nerve blood flow in 
glaucoma/' Semin. Ophthalmol, Vol.14, no. 3, pp. 164-170 (1999); A. Harris et al, 
"Simultaneous management of blood flow and IOP in glaucoma," Acta Ophthalmol Scand, 
Vol. 79, pp. 336-341 (2001); and S. S. Hayreh, "Factors influencing blood flow in the optic 
nerve head," J. Glaucoma, Vol. 6, pp. 412-425 (1997). Hypoxia of the retina and ONH is 
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believed to be a factor in the development of ocular vascular disorders such as diabetic 
retinopathy, arterial venous occlusion disease, and glaucoma. See, K. R. Denninghoff et al 9 
"Retinal imaging techniques in diabetes," Diabetes Technol Ther., Vol.2, pp. 111-113 
(2000); E. Stefansson et al. 9 "Oxygenation and vasodilation in relation to diabetic and other 
proliferative retinopathies," Ophthalmic Surg., Vol. 14, pp. 209-226 (1983); A. Yoneya et al 9 
"Retinal oxygen saturation levels in patients with central retinal vein occlusion," 
Ophthalmology, Vol. 109, pp.1521-1526 (2002); and E. Stefansson et al 9 "Optic nerve 
oxygen tension in pigs and the effect of carbonic anhydrase inhibitors," Invest Ophthalmol. 
Vis. ScL, Vol. 40, pp. 2756-2762 (1999). The ability to obtain relative measurements of 
oxygen saturation in the human ocular fundus could aid diagnosis and monitoring of these 
and other disorders. For example, measurement of changes in retinal and ONH oxygen 
saturation under controlled conditions could establish relationships between oxygen 
consumption, blood sugar levels, and vascular autoregulatory function in diabetic 
retinopathy. Assessment of oxygenation in the ONH may facilitate early detection of the 
onset of glaucoma, a disease in which timely diagnosis is crucial for effective treatment. 

[0005] Measurements of oxygen tension (PO2) in the ONH have been performed 

using O2- sensitive microelectrodes inserted into the eye. See, e.g., E. Stefansson et al. 9 
"Optic nerve oxygen tension in pigs and the effect of carbonic anhydrase inhibitors," Imest. 
Ophthalmol. Vis. Sci. 9 Vol.40, pp. 2756-2762 (1999). Although this technique is accurate 
and can determine p02 distribution in three dimensions, its invasive nature limits its use to 
animal models and precludes clinical applications. Another technique involving injection of a 
phosphorescent dye has been used to study p02 in the retinal and choroidal vessels, as well as 
the microvasculature of the ONH rim. See, e.g., S. Blumenroder et al 9 "The influence of 
intraocular pressure and systemic oxygen tension on the intravascular p02 of the pig retina as 
measured with phosphorescence imaging," Surv. Ophthalmol., Vol.42, pp. S118-S126 
(1997). However, use of the dye in humans has yet to be approved. 

[0006] Imaging techniques based on spectral changes of oxygenated hemoglobin 

(Hb02) and reduced hemoglobin (Hb) have been employed in humans to assess oxygen 
saturation in the ocular fundus, and in retinal artery/vein pairs. See Yoneya et al 9 (2002); and 
J. M. Beach et al. 9 "Oximetry of retinal vessels by dual-wavelength imaging: calibration and 
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influence of pigmentation, " Appl Physiol., Vol. 86, pp. 748-758 (1999). These methods 
have been based most often on recordings at several discrete wavelengths chosen for their 
relative sensitivity to changes in oxygen saturation. See, M. Crittin et al, "Hemoglobin 
oxygen saturation (So2) in the human ocular fundus measured by reflectance oximetry: 
preliminary data in retinal veins," Klin. Monatsbl. Augenheilkd, Vol. 291, pp. 289-291 
(2002); F. C. Delori, "Noninvasive technique for oximetry of blood in retinal vessels,"^/?/?/. 
Optics, Vol. 27, pp. 1 1 13-1 125 (1998); J. B. Hickam et al, "A study of retinal venous blood 
oxygen saturation in human subjects by photographic means," Circulation, Vol. 27, pp. 375- 
383 (1963); J. Hickam et al, "Studies of the retinal circulation in man: observations on vessel 
diameter, arteriovenous oxygen difference, and mean circulation time," Circulation, Vol. 33, 
pp. 302-316 (1966); and J. S. Tiedeman et al, "Retinal oxygen consumption during 
hyperglycemia in patients with diabetes without retinopathy," Ophthalmology, Vol. 105, pp. 
31-36(1998). 

[0007] Full spectral methods, employing a continuous range of wavelengths, have 

been used to record the reflectance profile versus wavelength from the ocular fundus. See, F. 
C. Delori, "Reflectometry measurements of the optic disc blood volume, " in Ocular Blood 
Flow in Glaucoma. Means, Methods and Measurements, G. N. Lambrou, E. L. Greve eds., 
Berkely, CA, Kugler and Ghedini, pp. 155-163 (1989); and F. C. Delori et al, "Spectral 
reflectance of the human ocular fundus," Appl Optics, Vol. 28, pp. 1061-1077 (1989). Full 
spectral imaging technique has also been employed to measure oxygen saturation in retinal 
arteries and veins under various conditions. See D. Schweitzer et al, "In vivo measurement 
of the oxygen saturation of retinal vessels in healthy volunteers," IEEE Trans Biomed Eng., 
Vol.46, pp. 1454-1465 (1999); and D. Schweitzer et al, "A new method for the 
measurement of oxygen saturation at the human ocular fundus," Int. Ophthalmol, Vol. 23, 
pp. 347-353. (2001). Oxygen saturation in the ocular fundus has been mapped using Fourier 
transform spectral imaging. See, Yoneya et al, (2002). The full spectral technique 
employed most often uses a high resolution imaging spectrograph to collect the spectral 
information from a band of tissue in a single spatial dimension. The method acquires data 
rapidly and is applicable for use in human subjects. See Schweitzer et al, (1999); and 
Schweitzer et al, (2001). 
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DISCLOSURE OF INVENTION 

[0008] We have discovered a new method to analyze continuous spectral curves to 

determine relative hemoglobin oxygen saturation. The method uses spectral curves collected 
from a continous range of wavelengths from about 530 nm to about 584 nm, including 
spectra from transmitted or reflected light. Using isosbestic points and curve areas, a relative 
saturation index was calculated. With this method, noninvasive, in vivo measurement of 
relative oxygen saturation was using light reflected from blood vessels in the eye. This 
method could also measure oxygen saturation from other blood vessels that reflect light 
sufficient to give a clear spectra from the blood hemoglobin, e.g., skin, tongue, or intestine. 
This method was used in connection with hyperspectral imaging to generate two-dimensional 
maps of tissues indicating relative hemoglobin oxygen saturation. In particular, this method 
was used to map and measure relative changes in hemoglobin oxygen saturation in primate 
retinal vessels and optic nerve head in response to controlled changes in inspired oxygen and 
intraocular pressure (IOP). Changes in blood oxygen saturation can be monitored with this 
method for early detection of disease, e.g., diabetic retinopathy or glaucoma. This method 
could also be used to monitor oxygen treatments for wounds or burns. 

Description of the Drawings 

[0009] Fig. 1 illustrates the hyperspectral imaging system in relation to the fundus 

camera. 

[00101 Fig. 2 illustrates an optical diagram of the retinal hyperspectral imager. 

[0011] Fig. 3 illustrates the organization of spatial (x,y) and spectral (X) information 

in acquired image frames (left) and after conversion to band-sequential images (right). 

[0012] Fig. 4 illustrates the reflectance spectra of saturated blood (Hb0 2 signature, 

bold curve) and desaturated blood (Hb signature, thin curve) from retinal recordings. 

[0013] Fig. 5 illustrates, on the left side, single-band images (570 nm) of the optic 

nerve head (ONH) and vessels from hyperspectral images obtained for oxygen breathing 
(top) and intraocular pressure (bottom) experiments, with labels corresponding to the vessel 
type (A, artery; and V, vein), and to nasal (N) and temporal (T) aspects of the optic nerve 
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head; and with lines indicating the optic nerve head areas, white lines for rim areas and black 
lines for cup areas. The right side shows flourescein angiograms (venous phase) used to 
confirm the vessel type (A, artery, and V, vein) for each experiment. 

[0014] Fig. 6 illustrates the spectral curves from various locations in the eye while 

breathing room air and pure oxygen. 

[0015] Fig. 7 illustrates the spectral curves from various locations in the eye at room 

pressure and at higher intraocular pressure. 

[0016] Fig. 8A illustrates partial signature saturation maps of the oxygen breathing 

experiment before (left, room air) and during pure oxygen breathing (right). The low to high 
oxygen saturation is indicated by the progression inside the blood vessel from blue (pale 
grey) to green to yellow (bright white) to red (almost black). 

[0017] Fig. 8B illustrates full signature maps of the oxygen breathing experiment 

before (left, room air) and during pure oxygen breathing (right). The low to high oxygen 
saturation is indicated by the progression inside the blood vessel from blue (pale grey) to 
green to yellow (bright white) to red (almost black). 

[0018] Fig. 9 illustrates the relative saturation indices (RSIs) from the oxygen 

breathing experiment (pure oxygen breathing started after the second data point, time = 0) 
from retinal vessels (large symbols: artery (filled diamonds) and vein (filled rectangles)) and 
from optic nerve head (ONH) regions (small symbols: temporal cup (open circles), nasal cup 
(open triangles), and average over ONH (x)). 

[0019] Fig. 10A illustrates full signature saturation maps of the intraocular pressure 

experiment with the top three maps under normal IOP conditions (15 mm Hg), and the 
bottom four maps under increased IOP (60 mm Hg), representing (left to right) 1, 2, 3, and 4 
minutes after onset of increased pressure. The low to high oxygen saturation is indicated by 
the progression inside the blood vessel from blue (pale grey) to green to yellow (bright white) 
to red (almost black). 

[0020] Fig. 10B illustrates full signature saturation maps of the intraocular pressure 

o experiment with a compressed scale showing (left to right) 1, 2, 3, and 4 minutes after the 
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onset of increased IOP. The low to high oxygen saturation is indicated by the progression 
inside the blood vessel from blue (pale grey) to green to yellow (bright white) to red (almost 
black). 

[0021] Fig. 1 1 illustrates the relative saturation indices (RSIs) from the intraocular 

pressure experiment (60 mm Hg pressure started after the third data point, time = 0) from 
retinal vessels (large symbols: artery (filled diamonds) and vein (filled rectangles)); and from 
optic nerve head (ONH) regions (small symbols: temporal cup (open circles), nasal cup (open 
triangles), and average over ONH (x)). 

MODES FOR CARRYING OUT THE INVENTION 

[0022] The optic nerve head (ONH) and overlying vessels in cynomolgus monkey 

eyes were imaged using a fundus camera attached to a hyperspectral imaging system. Images 
were acquired with inspiration of room air and pure oxygen, and at controlled intraocular 
pressures (IOP) of 15 mm Hg (normal) and 60 mm Hg (sustained for up to 5 minutes). 
Changes in relative blood oxygen saturation in the vessels and ONH were assessed from 
reflectance spectra. Saturation maps were derived from contributions of oxygenated and 
deoxygenated hemoglobin spectral' signatures extracted from hyperspectral images. The 
results obtained with hyperspectral imaging were compared with known experimental 
outcomes. 

[0023] Pure oxygen markedly increased oxygen saturation in veins; increases in 

arteries and the ONH were smaller. The results obtained with hyperspectral image analysis 
agreed with known changes in oxygen saturation from breathing experiments. Raising 
intraocular pressure (IOP) reduced saturation in all structures and resulted in profound 
desaturation of arteries. During sustained high IOP, a rebound in saturation was observed in 
the ONH. Spatial maps clearly showed the saturation changes in arteries, veins, and 
surrounding tissues. 

[0024] Using this same method, relative oxygen saturation in blood vessels from 

other areas of the body could be measured both in one and two dimensions. If measured in 
two dimensions, as in hyperspectral imaging as discussed in U.S. Patent No. 6,276,798, a 
two-dimensional map can indicate the relative oxygen saturation. Measurement of oxygen 
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saturation can be used as an indication of disease or a predictor of disease, e.g., diabetic 
retinopathy, glaucoma, ulcer, etc. It can also be used to assess the effectiveness of hyperbaric 
oxygen treatment of flesh wounds or burns by monitoring the changes in blood oxygen in the 
blood vessels near the wound or burn. Since the hemoglobin under normal conditions occurs 
inside the blood vessels, this method can also be used to enhance the detail of the vasculature 
in a tissue when using hyperspectral imaging to produce spectral maps. 

Example 1 

[0025] Materials and Methods 

[0026] Animals. The use of animals in this study was approved by the Louisiana 

State University Health Sciences Center Institutional Animal Care and Use Committee and 
conformed to current standards in the use of animals in ophthalmic and vision research. 

[0027] Two cynomolgus monkeys with normal eyes were used. The monkeys were 

anesthetized, and their eyes dilated. The initial ophthalmologic examination included 
fluorescein angiography, color and red-free fundus photography, and slit-lamp examination 
of the fundus. To measure oxygen saturation of the optic nerve head (ONH) and paired 
retinal vessels, a contact lens was placed on the cornea to prevent drying, and reflectance 
hyperspectral imaging measurements as described below were obtained in one eye of each 
monkey. During imaging, pure oxygen was administered to one monkey to directly control 
blood oxygen saturation, and intraocular pressure (IOP) was controlled in the other monkey 
using methods described below. 

[0028] Systemic Oxygen Saturation. An ear oximeter probe (Ohmeda 3700, 

Wallingford, Connecticut) was placed on the monkey's earlobe to measure systemic 
oxygenation. A tracheal tube was positioned at the trachea and connected to a small-animal 
breathing chamber (Quantiflex; MDS Matrx Co., New York City, New York). The oxygen 
chamber was supplied through a pressure regulator from an oxygen tank at a rate of 3 L/min 
at atmospheric pressure. This procedure brought the oximeter reading to 100% saturation. 
Hyperspectral images were obtained as described below while the monkey breathed room air 
and during inspiration of pure oxygen. 
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[0029] Intraocular pressure (IOP). To raise IOP, a 27-gauge needle was inserted into 

the anterior chamber of the eye under slit lamp examination. The needle was connected to a 
500-ml reservoir containing saline solution with 0.1 ml gentamicin (40 mg/ml), 0.O3 ml 
clindamycin (150 mg/ml), and 4 ml dexamethasone (4 mg/ml). IOP was raised by elevating 
the reservoir. IOP was monitored by means of a tonometer (Tonopen XL; Medtronic, 
Jacksonville, Florida). Imaging was performed at normal IOP (15 ± 2 mm Hg) and high IOP 
(60 ± 2 mm Hg), close to the pressure needed to stop vessel perfusion. High pressure was 
maintained for no more than four minutes, while recordings were made at one minute 
intervals. 

Example 2 

[0030] Hyperspectral Imaging System 

[0031] Fundus Camera. The retina was illuminated with the internal tungsten aiming 

light of a fundus viewing camera (TRC-50vt, Topcon, Japan), similar to the procedure 
described in U.S. Patents No. 5,919,132; and 6,276,798. Images were acquired using this 
camera with an ophthalmologic lens and a c-mount through the vertical path of the camera. 
Hyperspectral images were obtained through the vertical viewing port using an imaging 
spectrograph and digital camera, as described below. Fig. 1 shows the components and 
position of the hyperspectral imager on the fundus camera. The image normally formed at 
the film camera port. During hyperspectral imaging, the image is redirected upward by a 
mirror. The imaging system is translated over the camera port by a linear actuator mounted 
below the imaging spectrograph and charge-coupled device (CCD) camera. A vertical 
mounting facilitated image scanning by maintaining the center of gravity of the moving 

r 

components over the line of travel. A sleeve held the system at the proper height to sample 
the focused image. The entrance slit of the spectrograph was placed at a conjugated image 
plane of the eye fundus with the aid of the lens and c-mount. 

[0032] Hyperspectral Imaging. The hyperspectral images were obtained by 

translating an imaging spectrometer and charge-coupled device (CCD) camera (model VNIR 
100, Photon Industries Inc., Stennis Space Center, Mississippi) across the fundus image, as 
shown in Fig. 2 The area of interest on the retina is imaged with the fundus camera (FC). 
The dotted lines in Fig. 2 show the light collection path only. An intermediate image (IM) is 
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formed at the slit (S) of the imaging spectrograph (IS). The spectrograph is drawn above the 
image for clarity. The output spectrum is focused on the sensor of the CCD camera. As the 
spectrograph and the camera are translated along the Y axis, the spectrum from points on 
consecutive lines of the image is recorded in a series of frames. Motion was controlled to 
create a 1:1 aspect ratio between adjacent pixels in the X direction and lines in the Y 
direction. 

[0033] The spectrograph employed a prism-grating-prism (PGP) architecture with 2.5 

nm spectral resolution (25-um slit) and a range of 410 to 950 nrn. Images of the back of the 
eye were acquired using the 35° viewing mode of the fundus camera. The image from the 
vertical camera port was focused onto the entrance slit of the spectrograph. The output 
spectrum was in turn focused onto the CCD image sensor. This arrangement caused the 
spectrum of all points along a line in the fundus image to be recorded in a single CCD frame. 
Frames contained a maximum of 1024 points per line and 1024 points per spectrum. 

[0034] If the highest spatial or spectral resolution was not needed, greater light 

sensitivity could be obtained by binning CCD pixels. In the examples described below, two 
spatial and four spectral pixels were binned together to give spectral images containing 512 
spatial points and 256 spectral bands. This resulted in sufficient light sensitivity of individual 
picture elements and sufficient spatial resolution to enable us to monitor oxygen-dependent 
spectral changes in vessels. The second spatial dimension was obtained by translating the 
imaging system at constant velocity in the direction transverse to the orientation of the slit. 
The translation system comprised two mounts attached respectively to the fundus camera and 
the spectrograph, and a servo-controlled actuator that provided linear motion between these 
parts; relative motion of this system caused the slit to remain in focus with the fundus image 
throughout the scan. This component is termed the focal plane scanner (FPS). The number 
of rows obtained in each hyperspectral image was equal to the number of frames acquired as 
the system was translated. The velocity of motion and the interval between frames was 
carefully adjusted so that adjacent pixels and adjacent rows of the image had the same spatial 
interval. Typically 100 rows were obtained for this study. 

[0035] Fig. 3 shows the data structure of the recorded spectral images. Each frame 

holds the spatial (X) and spectral (X) axis for each line of the acquired hyperspectral image 
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(left image), with successive lines forming the Z-axis in the stack of frames. A "band- 
sequential" hyperspectral image is obtained by rotation of the stack of images, interchanging 
the Z and X axes. After rotation, each frame contains a two-dimensional spatial image (right 
image) at a distinct wavelength in which intact structures are recognizable. 

[0036] Extraction of Spectral Curves. Band-sequential image sets were saved from 

the image acquisition software (HyperVisual™; ITD, Stennis Space Center, Mississippi) in 
ENVI image processor format (ENVI, Research Systems, Boulder, Colorado). Images were 
corrected for dark values by subtracting an image obtained after blocking illumination. 
Spectral curves were obtained in ENVI by scanning the intensity profile along the Z-axis of 
selected image pixels within the optical nerve head (ONH) border, corresponding to artery, 
vein, and surrounding ONH. For spectral curves, a five-point moving average filter was 
applied to individual curves of each time point, and the smoothed data were then averaged to 
obtain final curves that represent the spectral signatures obtained before application of high 
oxygen, after application of high oxygen, and before high IOP. Time points during high IOP 
were not averaged. 

Example 3 

[0037] Mapping Relative Oxygen Saturation. 

[0038] Reference Spectra for High and Low Oxygenation States: Relative saturation 

was assessed from amplitudes of the hemoglobin spectral signatures that were contained in 
the reflectance spectra from retinal blood. As saturation was decreased from a high to a low 
value, spectral minima at 542 and 577 nm from oxygenated hemoglobin (Hb0 2 spectral 
signature) were converted to a single minimum at 555 nm from deoxyhemoglobin (Hb 
spectral signature). No changes occurred at wavelengths where Hb02 and Hb spectral curves 
crossed (isosbestic points). These spectral features from reflectance recordings at high and 
low saturation are shown in Fig. 4. Although the sloping baseline produces a slight blue-shift 
of spectral minima, only the areas under curves are used in this method. In Fig. 4, reflectance 
spectra of oxygen-saturated blood (Hb0 2 signature, bold curve) and oxygen-desaturated 
blood (Hb signature, thin curve) are shown from retinal recordings. The Hb0 2 curve contains 
two minima corresponding to wavelengths of peak light absorption. The Hb curve contains a 
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single broad minimum. The sloping baseline causes minima wavelengths to be shifted 
slightly to shorter wavelengths. Vertical lines extend from the axis to wavelengths where 
HbC>2 and Hb have equal reflectance and absorbance (isosbestic points of hemoglobin in 
distilled water). Dotted lines connect pairs of isosbestic points. Regions I, II, and III are 
defined between isosbestic points: Region I between 530 and 545 nm; Region II between 
545 and 570 nm; and Region III between 570 and 584 nm. In each region, the area between 
the spectral curve and the dotted line is denoted as al, all, and alii, respectively. AI, All and 
AIII represent the areas under the dotted lines to the baseline in regions I, II, and III, 
respectively. Saturation maps were determined from region II (partial signature) and from the 
combination of regions I, II, and III (full signature). 

[0039] Isosbestic points at 530, 545, 570, and 584 nm were selected from recorded 

spectra. As seen in Fig. 4, the curve of saturated blood passes above the line that connects the 
points at 545 and 570 nm (region II). The curve moves toward the line and passes below the 
line as the blood becomes more desaturated. This area between the spectral curve and line 
connecting points at 545 and 570 nm (designated "all" in Fig. 4) is largest for 100% 
saturation and decreases, eventually changing sign, as the blood becomes desaturated. 
Changes in the total reflectance from different recordings were compensated for by dividing 
this saturation-sensitive area by the area between the line connecting points at 545 and 570 
nm and the baseline (All in Fig. 4). This area is proportional to the intensity of reflected light 
in the recorded spectrum and is not affected by saturation changes. A partial signature map 
of relative oxygen saturation was found from the ratio of these saturation-dependent and 
saturation-independent areas in region II (ail/AH). A partial signature index could be 
calculated using this ratio. The term partial signature refers to the use of only the region 
(Region II) of the spectrum between the second pair of isosbestic points. 

[0040] A second method, producing a full signature saturation map, used three 

regions of the spectral curve as shown in Fig. 4. The first region (Region I) is defined as the 
region between 530 and 545 nm; the second region (Region II) is defined (as above) as the 
region between 545 and 570 nm; and the third region (Region III) is defined as the region 
between 570 and 584 nm (Fig. 4). Spectral maps were produced that included all three 
regions to determine if a significant reduction in noise and increase in sensitivity could be 
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obtained by all three regions. Areas between the curve and the line connecting the isosbestic 
points in regions I and III were negative at high saturation and moved toward zero and even 
slightly positive for low saturation. This second map, referred to as afiiH signature map, was 
found by subtracting areas I and EI from area II (after each area was compensated for total 
reflectance differences as described above for region II). The full signature map then was the 
ratio all/All minus the two ratios al/AI and alll/AIII. The use of the full signature gave a 
larger range of values for the same change in saturation and tended to average noise to a 
greater extent. For each type of map, values representing low to high saturation were color- 
coded as blue, green, yellow, and red. Because spectral changes were referenced to isosbestic 
points, this method minimized errors contributed by variation in the slope of the spectral 
baseline from different recording sites. 

[0041] Relative Saturation Indices. An index of the relative oxygen saturation (RSI) 

was determined from separate regions of the hyperspectral image containing artery, vein, and 
selected areas of the ONH (see Fig. 5). For each region an average spectral curve was 
determined, and then the RSI was calculated by the same method described above for 
individual spectral curves. 

Example 4 

[0042] Results of Either Pure Oxygen or Increased Intraocular pressure on 

Spectral Curves. 

[0043] Fig. 5 shows the area of the optical nerve head (ONH) obtained from the 570 

nm band in the hyperspectral image for the oxygen concentration images (top left) and the 
variable IOP images (bottom left). Confirmation of the vessel type was done by fluorescein 
angiography (images at right from the venous phase) for each experiment. Relative 
saturation indices in Figs. 9 and 1 1 were determined from retinal vessels marked as A (artery) 
and V (vein), and from the ONH inside areas bounded with white lines (rim) and black lines 
(cup) as labeled in Fig. 5. Nasal and temporal aspects of the ONH are labeled N and T 
respectively. 
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[0044] Spectral Signatures. Figs. 6 and 7 show a portion of the reflectance spectra 

between 450 - 600 nm containing the hemoglobin signature (spectral curves) from retinal 
artery, vein, and nasal and temporal ONH under various experimental conditions. Increased 
oxygen saturation is indicated in these plots when the experimental spectrum changes to more 
closely match the Hb0 2 signature of Fig. 4, with stronger minima at 542 nm and 577 nm. 
Desaturation is indicated when the curve more closely resembles the deoxyHb signature 
having a single spectral minimum. Higher reflectances at the longer wavelengths result 
mainly from weaker light absorption at these wavelengths by choroidal pigments in the 
fundus. The spectral curves of both Figs. 6 and 7 showed the expected Hb0 2 signature in 
arteries and the mixed Hb-HbC>2 signature in veins under room air conditions. Switching to 
pure 0 2 strengthened the Hb0 2 signature of both types of vessels, as seen in Fig. 6. If retinal 
arterial saturation is closely matched to the systemic saturation (95-97%), observed increases 
in the Hb0 2 signature in the artery represent only a 3-5% increase in saturation. Oxygen 
leakage from the ophthalmic artery could cause the retinal artery saturation to be lower than 
systemic levels. In that case, the response seen in the artery may represent up to an 8% 
increase in saturation. Since a fixed leakage rate would result in more or less arterial 
saturation depending on the flow rate, evaluation of retinal arterial saturation could 
effectively probe changes in blood flow at the major vessels supplying blood to the inner 
retina. The proportionately stronger Hb0 2 signature observed in veins corresponds to 
significantly larger increases in venous saturation. This effect was noted previously and was 
attributed to inhibition of the desaturation of capillary blood in the presence of high plasma 
p0 2 . See J. Hickam et al (1966); and J.M. Beach et al (1999). 

[0045] Oxygen breathing. The effect of inspired 0 2 concentration is shown in Fig. 6. 

The artery (top left) showed a small increase in the Hb0 2 signature with pure 0 2 , relative to 
room air. In the vein (top right), this increase was markedly larger. Inspiration of pure 0 2 
raised total reflectance, as shown by the greater spectral amplitude. In the nasal and temporal 
ONH (bottom left and right), pure 0 2 increased the Hb0 2 signatures, but not to the degree 
observed in the vessels. The larger increase was in the nasal ONH spectrum. Pure 0 2 also, 
increased total reflectance from the ONH. All spectra from the ONH showed an increased 
baseline slope because of higher reflectance at red wavelengths. As expected, the overall 
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results show increased oxygen saturation in both the large vessels and the ONH 
microcirculation with increased concentration of inspired O2. 

[0046] Intraocular pressure (IOP). The effect of increased IOP on oxygen saturation 

as seen in the spectral curves is shown in Fig. 7. In the artery (top left), high IOP sustained 
for 5 minutes gradually converted the Hb0 2 signature to an Hb signature. Top left in Fig. 7 
shows artery at normal IOP and at 1,3, and 5 minutes after IOP was increased to 60 mm Hg. 
Top right shows a vein at normal IOP and 1 minute after pressure was increased to 60 mm 
Hg. Bottom left indicates the nasal ONH. Bottom right indicates the temporal ONH. For the 
high IOP curves, the dotted line represents 1 minute and the solid line represented 5 minutes 
after IOP was increased to 60 mm Hg. The data indicate that a deep level of desaturation 
occurred in the artery. In the vein (top right), the normal IOP curve shows a weak Hb0 2 
signature. Within 1 minute after the onset of high IOP, however, the curve was converted to 
a strong Hb signature. These results suggest that high IOP causes desaturation of the retinal 
blood supply in both arteries and veins. Increased IOP resulted in only modest increases in 
total reflectance. 

[0047] In the ONH (Fig. 7, bottom left and right), Hb0 2 spectral signatures were 

present at low IOP. One minute after IOP was increased to 60 mm Hg, the amplitude of the 
signature decreased. At 5 minutes, the nasal ONH curve was nearly parallel to that at 1 
minute, whereas the temporal ONH curve showed some small restoration of the Hb0 2 
signature. These results show that high IOP reduced saturation in the ONH microcirculation 
but to lesser degree than in the retinal circulation, and suggest that saturation was partially 
restored in some regions. 

Example 5 

[0048] Effect of Oxygen Breathing and IOP on Relative Oxygen Saturation 

[0049] Responses to Oxygen Breathing. Figs. 8A and 8B show spatial changes in the 

relative saturation of ONH structures during room air breathing (left) and 2 minutes after 
switching to pure oxygen (right). Fig. 8A shows maps using the partial signature method, 
while Fig. 8B shows maps using the full signature method. The partial signature maps (Fig. 
8A) reveal saturation differences; however, structures such as the large vein are more clearly 
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delineated during high saturation in the full signature maps (Fig. 8B, right panel). In both 
figures, increasing haemoglobin oxygen saturation is indicated by the progression inside the 
blood vessel blue (pale grey) to green to yellow (bright white) to red (almost black). 
Temporal to nasal orientation in each map is top to bottom. These results show that better 
definition of the changes was revealed in the full signature maps. Accordingly, the full 
signature method was used to map IOP saturation changes below, and to determine the RSIs 
from vessel and ONH areas. 

[0050] Under room air conditions, high saturation areas included outlines of arteries 

out to the ONH boundary. These vessels continued outside the ONH with a different 
saturation code. During pure oxygen breathing, saturation increased in the arteries, and new 
areas of high saturation appeared where veins were located. The ONH tissue surrounding the 
vessels, particularly on the nasal side, showed smaller increases in saturation. These results 
agree with the spectral changes shown in Fig. 6. Table 1 displays RSIs averaged over two 
tune points for room air breathing, or over five time points for pure oxygen breathing. All 
structures showed significant increases (P< 0.05) in the RSI during pure oxygen breathing 
(Table 1); the increase in the veins was nearly twice (factor of 1.9) that found in the artery, 
whereas smaller increases in the ONH (averaged over the cup and rim) were approximately 
half (factor of 0.52) that of the artery. A slow decrease in the saturation over time occurred 
RSIs calculated from the vein and ONH, but not from the artery (Fig. 9). Fig. 9 indicates the 
relative saturation indices (RSIs) from retinal vessels and ONH for the oxygen breathing 
experiment. RSIs were determined from vessel segments inside the ONH and from rim and 
cup regions as denoted in Fig. 5. The following symbols are used in Fig. 9: vessel segments 
(large symbols): artery (filled diamonds), and vein (filled rectangles); and ONH regions 
(small symbols): temporal (open circles), nasal (open triangles), and average over ONH (x). 
Breathing pure oxygen began immediately after the second data point (time = 0). 
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[0051] Responses to High IOP. Hyperspectral imaging showed good repeatability, as 

is evident in the full signature saturation maps in Figs. 10A and 10B. Fig. 10A, top row, was 
from repeated recordings during low IOP (room pressure, 15 mm Hg). High saturation 
appears at artery locations within the border of the ONH and in the ONH tissue surrounding 
the vessels. Changes in saturation at 1 minute intervals after switching to high IOP are 
shown in Fig. 10A, bottom row. The bottom row (left to right) represents 1, 2, 3, and 4 
minutes after the onset of high IOP (60 mm Hg). Fig. 10B represents maps using a more 
sensitive scale (a scale that spans a smaller range) showing (left to right) at 1, 2, 3, and 4 
minutes after the onset of high IOP. Low to high saturation is indicated by the progression 
inside the blood vessel from blue (pale grey) to green to yellow (bright white) to red (almost 
black). For each map, temporal to nasal orientation is right to left. 

[0052] The high saturation of the arteries and most of the ONH disappeared after 1 

minute. A gradual return of saturation over the temporal ONH cup was observed from 2 
through 4 minutes after IOP elevation. 

[0053] Relative saturation indices are given in Fig. 11 and Table 2 for the IOP 

experiment. Table 2 displays RSIs averaged over three time periods for normal IOP, and 
over two time points for high IOP. High IOP resulted in significant reduction (P< 0.05) of 
the RSI. for each structure. After 3 minutes of high IOP, the values from artery and vein were 
not significantly different from one another. In the temporal ONH cup, the RSI decreased 
initially, but then recovered 24% of its original normal IOP value within the 4-minute high 
IOP period. This phenomenon was not observed in other areas of the ONH. Fig. 1 1 indicates 
relative saturation indices from retinal vessels and ONH for IOP experiment. The following 
symbols are used in Fig. 11: vessel segments (large symbols): artery (filled diamonds), and 
vein (filled rectangles); and ONH regions (small symbols): nasal cup (open triangles), 
temporal cup (open circles), nasal rim (filled triangles), temporal rim (filled circles), and 
average over ONH (x). High IOP begins immediately after third point (time = 0). 
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Table 2. Relative Saturation Indices from IOP Experiment 


Condition 


Artery 


Vein 


Averaged ONH 


Temporal cup 


Normal IOP* 


0.210 + 0.008 


0.139 ±0.005 


0.101 ±0.001 


0.081 ±0.001 


High IOP T 


0.030 ± 0.01 


0.029 ± 0.01 


0.041 ± 0.002 


0.054 ±0.001 


Difference 


0.180 + 0.018 


0.110 ±0.015 


0.060 ± 0.003 


0.027 ± 0.002 



Unpaired samples of equal variance. *Average over three time points at normal IOP. 
f Average over last two time points at high IOP. *A11 differences are significant (P< 0.05). 
Values are means ± standard deviations. 
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[0054] HbC>2 signatures were also obtained from areas between vessels within the border 

of the ONH. Since the light probe is in the green-red spectral range, these readings were 
interpreted to be signatures of blood carried by the microcirculation near the surface nerve fiber 
layer. It is also possible that some of this signal resulted from light first passing through surface 
vessels and then returning through the microcirculation of the surrounding tissue. Pure O2 
strengthened the Hb02 signature in the ONH, but to a lesser degree than that observed in the 
vein, as expected if this signature represents the averaged blood saturation in the 
microcirculation. These results are the first report of measurements of oxygen saturation 
changes in the ONH microcirculation using non-invasive reflectance imaging. 

[0055] Under pure O2 conditions, the ONH and vessel reflectance at the hemoglobin 

absorption wavelengths was consistently greater than under room air conditions. This effect may 
be the result of vasoconstriction under high O2 that reduces the luminal blood volume in the 
surface vessels and, correspondingly, the perfusion of the microcirculation. The features of the 
spectral profiles of vessels and tissue are thus in agreement with changes anticipated when the 
vascular supply of O2 is increased. 

[0056] Since metabolic changes associated with progression of retinal disorders 

presumably alter the oxygen utilization in the tissues, venous saturation maps should be a 
sensitive probe for disease states. Saturation maps determined by assessment of the Hb and 
Hb02 spectral signatures, in particular the relative contributions of the Hb and Hb02 spectral 
peaks between isosbestic points, were able to monitor the venous saturation increases in response 
to breathing pure O2. Previous work estimated these increases in the range 8-23%. If changes of 
similar size are present during the state of hypoxia, maps drawn as indicated above should be 
able to isolate hypoxic areas when the scale is moved to operate over the lower venous saturation 
range. Calibration for different saturation ranges would make the maps more sensitive in low 
and high saturation regions. 
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[0057] Raising IOP to 60 mm Hg had essentially the opposite effect on blood saturation. 

At this IOP, the perfusion pressure is very low. Arterial desaturation could have resulted from a 
slowing or stoppage of flow caused by collapse of the vessel under pressure, during which time 
oxygen diffused from the vessel. The more rapid appearance of the Hb signature in the veins 
was likely due to lower initial saturation of venous blood. An interesting feature of the high IOP 
response was partial recovery of saturation in the ONH microcirculation while the pressure 
remained high. Saturation recovery was seen near the cup of the ONH, which was temporal with 
respect to the origin of the vessels. The full signature map reduced noise enough to allow good 
visualization of this recovery. Since the high IOP effectively occluded the surface vessels, the 
source of oxygen is most probably from deeper levels of the circulation, which includes the 
retrolaminar layer. Increased reflectance during high IOP can be explained by low blood 
volume, since high IOP would partially occlude the major surface vessels and vessels feeding the 
outer ONH microcirculation, causing this area to blanch. 

[0058] Hyperspectral Imaging. These results demonstrate the ability of hyperspectral 

imaging to measure relative changes in oxygen saturation of blood vessels, e.g., retinal macro- 
and microcirculation. The usefulness of relative measurements of the oxygen saturation for 
assessing the vascular response to controlled changes in oxygen supply and utilization is evident 
from these data. 

[0059] The present hyperspectral imaging technique enables spectral quantitation to be 

carried out over two dimensions on the ONH, allowing regional changes in saturation to be 
identified. Different saturation color codes from retina outside the ONH were obtained. This 
difference may reflect disparate amounts of light being scattered into vessels from the pigment- 
free ONH and pigmented retina. In addition to the current method of curve integration, other 
spectral quantitation methods, such as curve fitting, can be employed. Significantly faster 
recording techniques would be better to achieve a more clinically acceptable method for 
mapping spectral information on the ocular fundus. 
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[0060] Hyperspectral imaging used with the current method can provide a much needed 

diagnostic tool for prevention and treatment of retinal disorders. The desired goal is the 
successful application of therapeutic interventions before irreversible damage occurs. One 
potential gain for detecting abnormalities in the oxygen saturation response is significantly 
earlier diagnosis of glaucoma. It is presently believed that autoregulation of blood flow is 
impaired in glaucomatous disease, possibly as a result of anatomical vascular impairment of the 
retina and the ONH. With this technique, problems in autoregulation could be diagnosed at an 
early stage, during the pre-onset stages of early phase glaucoma. In addition changes in oxygen 
saturation caused by other problems, e.g., diabetic retinopathy, hypertension, sickle cell anemia, 
and vascular diseases, can be detected by this method. In addition, this technique can be used to 
monitor oxygen saturation changes or blood flow in blood vessels from other body tissues, e.g., 
the skin, tongue, or intestine. The technique can also be used to assess skin disorders that might 
affect blood flow, e.g., a wound, a burn, or rosacea. The technique can be used to identify and 
locate major blood vessels in various regions of the body. 

[0061] The complete disclosure of all references cited in this specification are hereby 

incorporated by reference. Also incorporated by reference are the complete disclosures of the 
following: B. Khoobehi et al. y "Hyperspectral Imaging of oxygen saturation in the optic nerve 
head, retina, and choriod," Abstract presented May 7, 2003 at Association for Research in Vision 
and Ophthalmology; B Khoobehi et al., "Non-invasive measurement of oxygen saturation in 
optic nerve head tissue," Proc. SPIE, vol. 5325, pp. 104-110, Optical Diagnostics and Sensing 
IV; June 2004; and B. Khoobehi et aL, "Hyperspectral imaging for measurement of oxygen 
saturation in the optic nerve head," Investigative Ophthalmology and Visual Science, vol. 45, pp. 
1464-72 (2004). In the event of an otherwise irreconcilable conflict, however, the present 
specification shall control. 
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